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Abstract 

We have investigated the zero-temperature equation of state (EoS) for infinite nuclear matter within 
the {a — cj) model at all densities tib and different proton-neutron asymmetry rj = (N — Z)/{N + Z). 
We have presented an analytical expression for the compression modulus, and found that nuclear 
matter ceases to saturate at rj slightly larger than 0.8. Afterward, we have developed an analytical 
method to determine the strong coupling constants from the EoS for isospin symmetric nuclear 
matter, which allow us to reproduce all the saturation properties with high accuracy. For various 
values of the nucleon effective mass and the compression modulus, we have found that the quartic 
self-coupling constant G4 is negative, or positive and very large. Furthermore, we have demonstrated 
that it is possible (a) to investigate the EoS in terms of tib and r;; and (b) to reproduce all the known 
saturation properties without G4. We have thus concluded that the latter is not necessary in the 
{a — iu) model. 

PACS numbers: 

Kejrwords: Quantum Hadrodynamics, Relativistic Mean-Field Theory, Nuclear Forces, a lu p mesons. Asym- 
metry Energy. 
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I. INTRODUCTION 

Since the beginning of the 20th century, nuclear physics 
has dealt with the physical properties, structures and re- 
action mechanisms of the nuclear matter. It is a relatively 
recent branch of physics, which was developed intensely 
after the discovery of the atomic nucleus. One of the rea- 
sons for this is that its Coulomb field is shielded by the 
field of the bound electrons in the atom and therefore, 
it was difficult to prove the existence of the atomic nu- 
cleus in our environment. Furthermore, the dimensions 
of an atomic nucleus are so small that one could describe 
it technically and from there quantitatively only for the 
first time in 1911 using Rutherford's scattering experi- 
ments. This signalled the birth of nuclear physics. 

Nowadays, modern scattering experiments provide di- 
rect evidence for the following empirical properties of nu- 
clear matter: (a) The masses of the proton and neutron 
in the nuclear medium are approximately the same. The 
force between two protons is equal to the force between 
two neutrons (i.e., charge independence of the nuclear 
force) , but it depends strongly on the relative spin orien- 
tation. These properties can be proven by studying the 
inelastic scattering of protons and neutrons by deformed 
nuclei or by studying the mirror nucleus. Therefore, one 
introduces the concept of isospin symmetry by consid- 
ering the proton and the neutron as two different quan- 
tum states of the same particle, the nucleon. (b) The 
nuclear force has a very short range (approximately tt- 
meson Compton wavelength). This property provides a 



good approximation of the nuclear radius. The estimated 
value by the scattering of low-energy electrons from nu- 
clei is Q, i 



To ~ 1.16 fm. 



(1) 



'Electronic address; |Anis.Dadi@mpi-hd.mpg.de| 



which is 10^ times smaller than the atomic radius, (c) 
There is no attractive center, which produces the (strong) 
binding potential. The exact form of the nuclear po- 
tential is not precisely known at present, while in the 
atom, the Coulomb potential (oc l/r) is responsible for 
binding the electrons to the nucleus, (d) A nucleon may 
only interact (strongly) with a fixed number of its near- 
est neighbors and next-nearest neighbors. If one adds a 
further nucleon, only the nuclear volume becomes larger, 
not the binding energy per nucleon. From this it fol- 
lows that nuclear matter is a saturated system, (e) The 
medium is homogeneous and isotropic; that is, all direc- 
tions are equivalent and there is no direction preferred to 
others. This rotational symmetry is, as see in what fol- 
lows, very helpful for investigating the equation of state 
(EoS) for nuclear matter and generally in theoretical nu- 
clear physics, (f) The mean free path of the nucleon is 
too much larger than the nucleon size or even larger than 
the nuclear medium itself, so that one can apply a model 
of "independent particles". 

From these empirical properties, the nuclear medium 
can be compared with a liquid drop, where nucleons 
play the role of "molecules" in normal liquid. This phe- 
nomenological ansatz is known as "liquid-drop model" , 
whereby the behavior of the total binding energy B is 
determined as function of the total number A of nucleons 
and its different value for different isotopes. For a given 
nuclear medium, such as an atomic nucleus of volume V 
and total mass M, consisting of Z bound protons and 



N ~ A ~ Z bound neutrons of free masses nip and rrin, 
one assumes B to be approximately the sum of different 
types of energies, 

-^ -^volume 1 -^asymmetry i -^surface i -^Coulomb i ■■•; 

a -A-a .^Jl^^-a ^-A^'- 
ac-Z{Z-l)A-^/^ + ..., (2) 

which is equal to the mass defect {M — Zirip — Nm^)c^ . 
This is known as the semi-empirical mass formula Q- 
The coefficients Ov, flsym, ^surf, and ac are constants 
that can be calculated empirically by fitting to exper- 
imentally measured masses of different nuclei. In this 
work we focus only on the first and second terms. The 
first corresponds to the volume energy, which results from 
the mutual attraction between nucleons. Because the at- 
traction is short-range and rotationally symmetric, the 
binding energy per nucleon E = B /A is considered to be 
independent of A, and thus B is proportional to A. So 
the first term gives an explanation of the experimentally 
observed constant density of nucleons ne, whereby the 
volume V is proportional to the number A = V -n^ oi nu- 
cleons. If we simply estimate V as 47rrQ/3, using Eq.(IT]), 
the saturation density in the nuclear interior amounts an 
extremely high value of 



{n\ 



sat 
P 



') « 0.153 Nucleon/fm^ 



(3) 



where rip = Z/V and Un = N/V are the densities of 
protons and neutrons. Furthermore, the empirical value 
of the binding energy per nucleon is known with good 
accuracy from experimental data [l|, |^ : 



E. 
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-16.3 MeV. 



(4) 



This value can be obtained from the semi-empirical mass 
formula ([2]) by considering an isospin symmetric nu- 
clear matter («„ = Up) in the thermodynamic limit 
{A,V — )► oo), so that all terms vanish except the vol- 
ume energy. 

The second term in expression ([2|) has a quantum me- 
chanical meaning, namely, the nuclear medium can be 
considered as a quantum degenerate Fermi gas, in which 
all energy levels are occupied separately by protons and 
neutrons up to the Fermi energy. From the Pauli exclu- 
sion principle, each level allows just two particles of the 
same type with opposite intrinsic spin orientations. Ow- 
ing to this prevention, the more nucleons are added, the 
higher are the energy levels that they have to occupy, in- 
creasing the total energy of the medium and decreasing 
the binding energy (the latter will be maximum when 
protons and neutrons occupy the lowest possible levels). 
If there is an unequal number of protons and neutrons 
(usually there are more neutrons than protons), the en- 
ergies of the highest occupied proton and neutron levels, 
which are the Fermi energies, are unequal: 



(pV 



2?7lr 



^ 



2m„ 



In that case, a contribution to the binding energy should 
assure that protons and neutrons have the same Fermi 
energy; otherwise, the energy difference would be bal- 
anced through /3 decay of neutrons into protons or vice 
versa. This contribution is the asymmetry energy. Thus, 
the asymmetry energy provides for an equilibrium be- 
tween proton and neutron number, vanishes for n^ = ?^p, 
and lowers the total binding energy B by increasing the 
difference Un ~ Up . The empirical value of the asymme- 
try energy coefficient at the saturation density ^ is also 
well known H: 



<^^ « 32.5 MeV. 



(5) 



The next important quantity in the study of nuclear mat- 
ter properties is the nucleon effective mass m^ , which re- 
sults from the mass defect of the nucleon owing to the in- 
teraction with its nearest neighbors. From experiments, 
one determines m^ by measuring the density of states 
of the nucleon and the single-particle energy levels in nu- 
clei |5|. However, a precise value is not known at present. 
So far as known, the range at saturation density, far away 
from the Fermi surface and for n,i = rip, is 



II 



0.7 <m^/mN <0.8. 



(6) 



The compression modulus /C is also an important quan- 
tity in the study of nuclear properties, neutron star, su- 
pernova collapse and heavy-ion collisions. It is connected 
with the vibration process in the nucleus, which can be 
excited to oscillate through inelastic scattering of a par- 
ticles or through absorption of 7 quanta. Thereby, there 
are radial and bending vibrations, in which the deforma- 
tion of the nucleus changes with the oscillation periods. 
Concretely, /C can be extrapolated from experimental 
data on the strength function distribution of the isoscalar 
giant monopole resonance (GMR) and the isoscalargiant 
dipolc resonance (GDR) in different heavy nuclei |8l-[l0|. 
Like the nucleon effective mass, a precise value of the 
compression modulus is not known at present, but in re- 
cent years it has been estimated to be in the range pi)| - KL3| 

230.0 < /C < 250.0 MeV. (7) 

Furthermore, the value K. ~ 234.0 MeV, which we sup- 
port for the rest of this report, has been estimated by 
Myers and Swiatecki [2| using Thomas-Fermi model cal- 
culations. Youngblood et al. have reported [Ol that 
K, = 231 ± 5 MeV by measuring the strength function 
distribution of the GMR in ^OZr,"^ Sn,"4 Sm, and ^ospb 
using inelastic scattering of 240 MeV a particles at ex- 
tremely forward angles. 

The aim of theoretical nuclear physics is now to de- 
velop mathematical models, that provide an accurate de- 
scription of the observed physical properties in nuclear 
medium. On one hand, one tries to reproduce the known 
empirical data as optimally as possible, and on the other 
hand, one tries to make qualitative and quantitative pre- 
dictions on the physical behavior of the nuclear medium, 
which are unknown so far. Because the formulation of 



the quantum mechanics and the relativity, different the- 
oretical approaches have been developed, which can be 
divided into two different models: 

i. Microscopic models based on realistic nucleon- 
nucleon interaction. 

ii. Effective models based on a phenomenological 
parametrization of the interaction potentials. 

Numerical solution of realistic models, such as the (Dirac- 
) Briickner-Hartree-Fock theory, is very complicated, so 
they cannot be applied to finite nuclei. Rather simple and 
successful approaches generally offer the phenomenologi- 
cal models, such as the Skyrme model, the effective global 
color model of QCD, the Nambu-Jona-Lasinio model, 
and the relativistic {a — oj) model. The latter is the sub- 
ject matter of the present report, whose standard ver- 
sion is known as Walecka model [14| or quantum hadro- 
dynamics (QHD). This effective model provides an ac- 
curate description of the nucleon-nucleon interaction at 
large distances but breaks down at short distances, be- 
cause it treats baryons and mesons as elementary degrees 
of freedom and, of course, at short distances quarks and 
gluons manifest themselves. In contrast, the theory of el- 
ementary degrees of freedom of quarks and gluons, that 
is, QCD, is simple at short distances but complicated 
at large distances (in the confinement region). The pur- 
pose that allows a transition from QCD out into QHD is 
known as "hadronization," which is, however, not solved 
so far. 

The present article is organized as follows: In Sec|lll we 
calculate the EoS for a static and infinite nuclear matter 
at finite temperature. Afterward, we focus on the zero- 
temperature EoS by investigating the case of asymmetric 
(rip 7^ Tin), isospin symmetric nuclear (n^ = Hp) and 
pure neutron matter (rip = 0) at all densities. In Sec lIIIl 
we show how to determine the strong coupling constants 
analytically from the EoS, which allow us to reproduce 
the saturation properties of isospin symmetric nuclear 
matter discussed earlier. That is 

Input parameters : Saturation properties. 
Output parameters : Coupling constants, 

which is exactly the opposite of what one usually does. 
Our conclusions and summary follow in Sec. IIVI We 
try to use precise values of physical quantities from the 
recent 2010 CODATA [H]; this would be important to 
minimize the standard deviation errors of the calculated 
coupling constants from the empirical values and vice 
versa. Throughout we use the natural Heaviside-Lorentz 
units 

h = c = kB ~ ^, 

where h is the reduced Planck's constant, c is the velocity 
of light and ks is the Boltzmann constant. The product 
of h and c has the dimension of energy times length and 
the value is [15[ 

h-c= 197.326950073242 MeV • fm = 1. 



The value of the averaged nucleon free mass used is [l5| 
TON = (mp + ?7in)/2 = 938.91867973 MeV. (8) 

II. THE (a ij) MODEL 

At the saturation density ([3]), despite the low value of 
the binding energy Q compared with the free mass of 
the nucleon ([8]), a relativistic description is necessarily 
attributable to the delicate cancellation between short- 
range repulsion and medium-range attraction. It can 
also be justified by comparing the Fermi momentum of 
the nucleon with its free mass, kp w 0.276?7i,n. The 
(ct — uj) model has been developed to reproduce the sat- 
uration properties of equilibrium nuclear matter and to 
describe the relativistic effects at higher densities such as 
in neutron stars [16[ . The exchange of w and a mesons 
between nucleons is known to provide the short-range re- 
pulsion and intermediate-range attraction in the nucleon- 
nucleon potential. The effective interaction is character- 
ized by the meson parameters such as their masses and 
coupling constants, which are "adjusted" to reproduce 
the saturation properties. This interaction can be ex- 
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Figure 1: (Color online) Effective nuclear potential in the {a — 
Lj) model illustrating short-range repulsion and intermediate- 
range attraction. 

pressed phenomenologically for equal numbers of protons 
and neutrons in terms of two Yukawa potentials. 



V,.^U\x\) = V^i\x\) + V^i\x\) 



2 p-m^\x\ 

Att \x\ 



47r 



-m„\x\ 



, (9) 



where go-, 5^, ma- and m,^ are the coupling constants 
and masses of the a and a; mesons, respectively, and \x\ 
is the relative distance between two nucleons. The mi- 
nus sign tells us that the a meson potential is attractive. 
Moreover, the heavier (lighter) the exchanged meson is. 



the more rapidly (slowly) decreases the potential with in- 
creasing distance \x\. Therefore, the range of the result- 
ing interaction Q is mostly determined by the masses of 
the exchanged a and a; mesons and the strength by their 
coupling constants g^ and g^j as shown in Fig. [TJ 



A. The {a — ui) Model at Finite Temperature 

We start our theoretical investigation of nuclear mat- 
ter properties by defining our system and then construct- 
ing the corresponding lagrangian density. We consider a 
static and uniform system of A bound nucleons (Z pro- 
tons and N neutrons) moving in a volume V, which is 
assumed to be very large so that we can neglect the sur- 
face effects. The total mass and the total binding en- 
ergy of the system are denoted by M and B, respec- 
tively. We denote its energy density by £ = M/V 
and its binding energy per nucleon hy E = B/A = 
S/nB — "T-N, where tt-b denotes the density of nucleons, 
Ub = {Z + N)/V = rip + Tin- Each nucleon field depends 
on space x = {xi, X2, xa) and time t = x^ and represented 
by a Dirac-spinor ip{x,t) ~ (p, n)* as an isospin doublet 
state of proton and neutron, V' = (N t, N J,)*. 
The internal structure of nucleons, a and w mesons are 
not taken into account but described as relativistic bary- 
onic and mesonic degrees of freedom. The a meson is 
isoscalar-scalar and the lo meson is isoscalar- vector (i.e., 
both are neutral with total spin J = and 1, respec- 
tively). As a consequence, they do not have isospin- 
dependence in the lagrangian density, and therefore they 
do not contribute in the asymmetry energy. The asym- 
metry energy originates from the kinetic energy of nu- 
cleons (as discussed in the Introduction) and from their 
interactions itself (osym = awn + Oint)- To obtain the 
empirical value ([5]), we must introduce the p meson into 
the lagrangian density of the theory, which is an isovec- 
tor meson and couples to the isospin density (i.e., to the 
difference between the density of neutrons and protons, 
iT-ii — Up). Thus, the lagrangian density can be written as 



£(2)=^ 



7*" ( id^ - .9wXm - 2-9p'^^M ) - ('^n - g^ 



ij 






-mlx^.X^ 






(10) 



respectively. The superscript "(2)" above C is related 
to the quadratic mass term m'^cfp' of the (free) u meson 
and, as we see later, leading to a linear expression be- 
tween the scalar a potential and its scalar density. That 
is why £(^' is known as the lagrangian density of the 
"linear {a — lo) model" [1J|- In particular, we will see, 
that this linear version cannot reproduce correctly the 
nucleon effective mass ^ and the compression modulus 
([7]). Therefore, a nonlinear extension of £^^-' introducing 
scalar self-coupling terms (counterterms) is necessary to 
reproduce all the five empirical values [Eqs. (l3])-(l7|] and 
to describe the ultrarelativistic limit of the EoS in the 
high-density approximation (above 0.6 fm~'^). Let us 
then replace the mass term of the free a meson with the 
following nonlinear scalar potential. 



i^^<^^^f/(0) = %g!+(^+(^, (11) 



containing, in addition, a three- and four body interac- 
tion. The form of this potential is equivalent to that 
proposed by Boguta and Bodmer [131 where 



6= -;- and c= -;-. 



The lagrangian density pop will be then be written as 



£(^) = V^ 

{9a<j^? (.9.0)' {g. 



7^ ( id^ - g^XM - 29pTb^ ] - {m^ - g„c, 



^ 



2G. 



3G3 4G/+2^^^^^^+^^'^^^ 



2G, 



_9p_ 

2G ''^^ -^""^ 



5^6^ - -B^.B'"-. 



(12) 



Here we define the following five coupling constants 
as phenomenological input parameters of the nonlinear 
model 



Gcr 



"j G(^ 



!f , Gp = ^, Gg and G4, (13) 



Here r denotes the nucleon isospin matrix; cj} is the field 
of the scalar a meson; Xfj. and F/^^^ — d^^Xv ~ d^Xti. ai'S 
the field and strength tensor, respectively, of the vector 
w meson; 6^, B^^^, = dfj^b^, — d^bp^, nip and gp are the field, 
strength tensor, mass, and coupling constant, respec- 
tively, of the vector-isovector p meson. This lagrangian 
density is invariant under the global SU(2) isospin sym- 
metry and the f/(l) gauge symmetry. It contains the 
Yukawa couplings g^-iptjj'ip and guj'4'l'^Xfj.''P oi the nucleon 
field ip to the scalar and vector meson fields and x^. 



with dimG^ = dimG(^ = dim Gp = [Mass ^], dimGg — 

! 
[Mass^^] and G4 is dimensionless, because dim£(^^ ^ 
[Mass"*]. From the Hamilton's principle of least ac- 
tion (Euler-Lagrangc equations) with respect to ^ = 
{(l},XtJ.,bp,ip} 






(14) 



we derive the field equations of the hadronic fields, which 
are 



9a 



Go 

^2 



Ua ,2 I .Vcr ,3 



d^d>^ + ^]cf>+ ^<p' + ^4>' = ff.V'V', (15a) 



G: 



CjA 



Q^P^^U ^ ^^. ^ g^^^.^^ 









gpi^Tj"^, 



1 



7^ ( idf, - g^Xfi - -^gpTbf, ) - (ton - go 



(15b) 

(15c) 

^' = 0. 
(15d) 



We now introduce the path integral formalism into sta- 
tistical physics to calculate the EoS. First, the motion of 
the nucleon in the volume V is described by the proba- 
bility amplitude of finding it at the position x' at time t' , 
when we know that it was located at point x at time t. A 
path integral representation of the probability amplitude 
can be obtained by discretizing the time interval i' — i up 
into L infinitesimal segments of fixed duration At and 
then taking the continuum limit {L — >■ cx)) at the end of 
the calculation. Without getting into details, the result 



(V'',i>,i) = (V.|e-^^'" (*'-*) IV-) 



VipVipV(t)Vx^.Vh^ exp 



i dt (fxC^'^^ 



, (16) 



where H^^^ is the Hamiltonian of our hadronic system. 
Second, one of the substantial tools of statistical physics 
is the grand canonical partition function 



Z(r,V^,Ai)=Tre 



-li[H-t,A) 



(17) 



from which all thermodynamical quantities (energy, pres- 
sure, compression modulus, etc.) can be derived. Here 
A is the particle number operator of nucleons and // is 
the nuclear chemical potential. By considering the nu- 
clear medium as a quantum degenerate Fermi gas, the 
occupied levels can be seen as continuous, because the 
possible values of the quantized wave number fc are very 
close to each other (Afc -^ 0, analog to the electronic 
band structure). Therefore, one cannot only consider 
the trace in Eq. (|17p as a sum over all microstates, but 
rather as a phase-space integral; that is. 



PV(V'|e"^(^"^^)|^). (18) 



Z{T,V,ii) 



If we compare now Eq. (|T6|l with Eq. (|T8l) . we see that 
the probability amplitude ()16p is expressed in Minkowski 
space-time, which is not a normed space (i.e., its inner 
product is not positive-definite). In contrast, the par- 
tition function (fTSj) is defined in phase space F, which 
is analog to Hilbert space, and this latter is a normed 
space. So, to unify both expressions, we integrate the 
lagrangian density C^'^> along the imaginary time axis. 



t — > — ir. take r ~ 1/T as "inverse temperature," and 
make the substitution H — > H^) -\- ^A. The partition 
function (fT5|) can thus be expressed as a path integral in 
Euclidean space (.t^,x) = (r, Xi,af2, a^s) to give 



exp 




£(*) + /xV'tV } ■ (19) 



To evaluate Z and to solve the field equations (fT5)) . we 
use the relativistic mean-field (RMF) approximation un- 
der the assumption that the nuclear medium is static 
and uniformly extended. Notice that this approximation 
is not directly justified to simplify the calculations con- 
siderably, but rather could be seen as a consequence of 
the observed empirical properties of nuclear matter men- 
tioned in the Introduction. The RMF approximation ne- 
glects the fluctuations of the meson flelds and replaces 
them with their mean values, namely. 



a — y {(t) = ao = g^4)o, 
^t^ — > (^m) = '^o = 9uXo, 
Pf, — > {Pf,) = Pa= gpbQ, 



(20a) 
(20b) 
(20c) 



which means that protons and neutrons move indepen- 
dently in classical mean fields (Jjq, XOi and bo- As a first 
consequence, the field equations (|15p simplify consider- 
ably to give 



ctq 



^"-a3 = a(V;^>) 



Go- 2 I ^(^ 3 r^ /T i\ '^°^ /-I 



3 G4 

^0 = G'„(V'7"V') = GujUb = Go; [rip + n^) 
Po = 7;Gp{in°r3i') = -Gp (rip - n„) , 



«7m^^ - 7o ( Wo + -^T^po 



(mN - (To) 



(21a) 
(21b) 
(21c) 

ij{x,t) = 0. 
(21d) 



In Eq. ()21ap we defined a Lorentz scalar density n^ of 
the nucleons, which depends on the Fermi momenta of 
protons and neutrons, k^ and fc^. In Eq. (|21cp we have 
T3 = +1 for protons and T3 = — 1 for neutrons. Further, 
the partition function (fT9|) simply yields 



Z ((To, Wo, Po) = const X 



exp 



Pi 



;3yjii^ + ^:u in lu 11^ 



2G, 



ZGrp 2Ct(7 "^^3 4Gr4 



dV 



k^dk 



In 1 



-I3{E, 



'o-f^oA 



-PiK+t^l'. 



(22) 



where we define, respectively, 



E*'':^'Jk'+m*^\ 



* dcf 



mN - o-Q, 



* dcf 



-Jzpo, 



(23) 
(24) 
(25) 



as an effective dispersion relation, effective mass, and 
effective chemical potential (or effective Fermi energy, 
E*p) of the nucleon. The d in Eq. ((22|) denotes the de- 
generacy factor, having the value 4 for isospin symmetric 
nuclear matter (7t,„ = rip) and 2 for pure neutron matter 
(rip = 0). Rigorously speaking, the partition function 
(|22|) must be expressed in terms of product over fc^ and 
k'p^ and then taking d — 2. From Eq. ([25|l it becomes 
clear that the isospin forces arising from the p meson 
exchange lead to add the quantity po/2 to the neutron 
Fermi energy and to substract the same quantity from 
the proton Fermi energy, while the repulsive forces 
arising from the w meson exchange lead to decrease 
the Fermi energy of the proton from {pP — po/2) to 
[pP — po/2) — Wo and of the neutron from (/i" + po/2) to 
(pP + po/2) — Wo- Moreover, the attractive forces arising 
from the a meson exchange lead to decrease the nucleon 
mass from ton to vri^ and this mass defect becomes 
more important by introducing the three- and four-body 
interactions (|lip . The scalar mean fields cro, i^o, and 
Po have the dimensions of mass; the mean effective 
quantities rrv^ and /ig are scalar in Dirac space and 
depend on the Fermi momenta of protons and neutrons. 

Now we are able to derive the EoS for infinite nuclear 
matter and to investigate its thermodynamics at given 
values of N and Z . In thermodynamic equilibrium, the 
meson-exchange process is reduced substantially, leading 
to a maximum entropy S, which is 



/ig, whereas the second one has the same effective mass 
but opposite chemical potential, —p^, so that the total 
proton number remains constant (from the beginning, we 
do not consider for example the /? decay). Further, the 
sum of both distributions is divided by the Lorentz factor 
Eq/iti^, which is an effect of the Lorentz contraction; 
the higher the density of protons is, the more important 
this relativistic effect will be. These interpretations are 
analogous to that for the scalar density of neutrons in 
Eq. (M . 

Similarly, the partial derivative of S with respect to Wq 
and po yields the "net" baryon densities of protons and 
neutrons, 



np(fcP,/.P,r) 



(27c) 



1 



1 /-^F 
/ k dk 

TT^ Jo Ve'^(^°-^5) + 1 e'^(^o+''o) + 1 



nn{kp,p'^,T) = 
At I "^ k'^dk ■ 



(27d) 



1 



,,3(B*-mS) 



1 6' 



fiiE^+f^o) 



B. Equation of State at Zero Temperature 

At zero temperature, the thermal Fermi distributions 
of baryons in Eas. (|27p give one, and of antibaryons give 
zero, because we still have k ^ kp ^ Eq ^ Ep = /ig. 
Thus, Eqs. (|27cp and (|27d|) yield the well known expres- 
sions for the baryon density. 



riB 



1 
3^ 
and the isospin density. 



{k^pf 



(Kf 



(^f)' 



(fc^)^ 



(28a) 



(28b) 



dS _ dS _ dS ^ 
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with S — lnZ((Tg, Wo, po)- So, using Eq. ([22]) . making 
the derivative of S with respect to uo and identifying the 
result with Eq. (|21a|) . we obtain the scalar densities of 
protons and neutrons at a given temperature T 
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As an application, for isospin symmetric nuclear matter 
at saturation density ([3]), we obtain 



f£ p — ' fS p 



259.147 MeV« 1.313 fm" 



(29) 

Now using Eqs. ([M]), dH^), and (H?]), we deduce the first 
zero-temperature EoS at given values of proton and neu- 
tron densities, namely, the nucleon effective mass 
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(30) 



The pressure can be deduced from the grand canonical 
potential J = -T ■ In Z(T, V,p) = -P -V to give 



with nP -I- n^ = Ua- [Eq. (|21ap ]. On the right-hand side 
of Eq. (|27ap . we see that the scalar density of protons is 
expressed as a sum of two thermal distribution functions 
for a relativistic Fermi gas of protons and antiprotons. 
The first has an effective mass mjj and chemical potential 
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(31) 



Further, the energy density E ~ M/V can be deduced 
from the Gibbs-Duhem equation, which is simply given 

by 



f (tib) = fi-HB - P{nB) 



(32) 



at zero temperature. This, together with Eqs. ()21bp and 
([25|l yields 



matter {Z ^ TV). So, if we write the binding energy per 
nueleon at given Z and N as 

E{nB, ri) == EisoiriB) + if ■ agym + 0{7f), 



then we get 



E{nB,l)-E{nB,0). 



(37) 



£ {np, Un) = -G„ {up + n^f + -Gp [n^ ~ Upf + -;^ + 
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Notice that all integrals over the Fermi momenta appear- 
ing in Eqs. ([50)1 . ([5T|) . and ([55)1 can be done analytically, 
and then they can be expressed solely in terms of fc^ and 
fc^ or, according to Eq. ([^5]) . the baryon densities Tip and 
rin- Further, from the semi-empirical mass formula ([2]), 
we see that the asymmetry energy is given by 






riB fixed 
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(34) 



where 77 = Up/uB = (n-n — np)/{na + rip) denotes the 
proton-neutron fraction. To calculate asym within the 
{a — Lo) model, it would be easier to express rip and n^ 
in term of 77 and ub , 



np = — (l-r?), 
'^n = ^(l + J?), 



(35a) 
(35b) 



and then make the second derivative of the energy density 
P3p with respect to 77. After some calculation, we obtain 



ki 



6EI 



+ -^"-B = flkin + flint- 



(36) 



The first term Okin arises from the difference between the 
Fermi energies of protons and neutrons (see Introduc- 
tion). The second a-mt arises from the isospin coupling 
between the p mesons and the nucleons. Both together, 
flkin + Ointj assurc that neutrons and protons have the 
same Fermi energy. 

Alternatively, we can easily prove from Eq. ^ , that agym 
is the difference between the binding energy per nueleon 
E(np, n,i) of pure neutron matter and that of isospin sym- 
metric nuclear matter. Namely, from Eq. (j35p we have 
r] = 1 for pure neutron matter, 77 = for isospin symmet- 
ric nuclear matter and < 77 < 1 for asymmetric nuclear 



We can thus compare the numerical result of Eq. (|57|) 
with that of Eq. ([M)l . In addition, we can investigate the 
evolution of the energy density ([55)1 at all values of the 
proton- neutron fraction 77 between and 1, which now 
depends solely on this parameter and tib. 

We close the theoretical framework by calculating the 
compression modulus. By starting from its definition and 
using the first law of thermodynamics, dB = —PdV, we 
obtain 
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(38) 



,T=0 



This simply means that the compression modulus is the 
slope of the pressure pip at saturation density and multi- 
plied by 9. Furthermore, /C may be expressed analytically 
in terms of 77 and the saturation density by evaluating the 
derivative of the scalar mean field cto in Eq. (|21aP with 
respect to 77b, or equivalently to kp, which yields the 
following expression in terms of 7/ and ub '■ 
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with 
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fcP = (37r277p)i/3, 

This, together with Eq. (|38p and the expression of the 
pressure pip yields K. for cold nuclear matter in terms 
of 77 and 7ig"* 
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(40) 



We point out that all quantities appearing in Eq. (|^ . as 
well as the derivative of the scalar mean field ([Mjl . should 
be evaluated at tib = nff*. 



Gcr, G^j, Gp, G3, and G4 used for the RMF4 model 
are number 6 of Table IIIII 



Results 



Before discussing the results, we should make clear the 
meaning of the following notations: 

• RMF2 means the relativistic mean-field approxi- 
mation to the linear {a — lo) model, that is, with- 
out three- and four-body self-interactions. The la- 
grangian density used for this is given by Eq. (|10p . 
There arc thus three coupling constants: G^r, Gi^, 
and Gp. 

• RMF3 means the relativistic inean-field approxi- 
mation to the nonlinear (a — cj) model using the 
lagrangian density (fT2|) but without the four-body 
self-interaction. The coupling constants are four: 
Go-; Gtj, Gp, and G3. 

• RMF4 means the relativistic mean-field approxima- 
tion to the nonlinear (a — uj) model by introducing 
the three- and four-body self-interactions. The cou- 
pling constants are five: G^, G^^, Gp, G3, and G4. 

One of our aims is to reproduce the five saturation prop- 
erties of isospin symmetric nuclear matter for each model 
by determining separately the three, four, and five cou- 
pling constants of RMF2, RMF3, and RMF4, respec- 
tively. Notice that only Eq. (|30p must be solved numeri- 
cally, that is. 



Figure [5] shows the nucleon effective mass ([5D|) as a 
function of the nuclear density only for isospin symmet- 
ric nuclear matter and pure neutron matter. At low den- 
sity, n^,kp — >■ 0, rri^ has almost the same behavior in 
the three models, because it converges to ton and then 
the nuclear medium has more tendency to an ideal gas 
behavior. Afterward, mj, decreases almost linearly with 
the density. At higher density, above 0.6 fm~'^, it con- 
verges slower to zero in both RMF3 and RMF4 than 
in the RMF2 model; the first-order Taylor expansion of 
Eq. (|30|) around m^/m-^ yields 



mN 



(41) 



G3G4 + Gam^{GA + G^m^) 



G3G4 + GomN(2G4 + 3G3mN) + GaGiG^k'p/n^ 



■m-li = f{nB,mli), 



and the other EoS dSTj), dSSj), §6^, and (gO]) can be 
deduced analytically. For example, using a Fortran- 
95 code for the secant method [2l|, the error we ob- 
tained is roughly 10~^ comparing with an exact solution. 
Throughout the rest of this section values of the coupling 
constants 

• Ga , Gu: , and Gp used for the RMF2 model are num- 
ber 5 of Table El 

• Go-, Gi^, Gp, and G3 used for the RMF3 model are 
number 15 of Tabic ITVl and 



which gives in the RMF4 model m^j w 238.3 MeV at 



ub 



1.0 fm and n„ 



N 

Further, the behavior 



of TTiJj for isospin symmetric nuclear and pure neutron 
matter is roughly the same at all densities owing to the 
approximative equality of the proton and neutron free 
masses and to the charge independence of the attractive 
force arising from the a meson exchange. At the satu- 
ration density (I3|), the coupling constants used here lead 
to m^ ~ 0.540mN in RMF2, 0.798mN in RMF3, and 
0.780mN in the RMF4 model. 
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Figure 2: (Color online) The nucleon effective mass as a func- 
tion of the nuclear density in the RMF approximation to the 
nonlinear {a — uj) models (RMF3 and RMF4) and the linear 
model (RMF2) at zero temperature. The solid lines corre- 
spond to isospin symmetric nuclear matter, whereas dashed 
lines correspond to pure neutron matter. 



Figure 3: (Color online) The binding energy per nucleon as 
a function of the nuclear density in the RMF approximation 
to the nonlinear {a - uj) models (RMF3 and RMF4) and the 
linear model (RMF2) at zero temperature. The solid lines 
correspond to isospin symmetric nuclear matter, dashed lines 
to pure neutron matter, and dot-dashed lines to the asym- 
metry energy. The empirical values ng"*, E^^^ and a|y^ are 
reproduced with an standard deviation error less than 10~*%. 
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The binding energy per nucleon is illustrated in Fig. [3] 
for both isospin symmetric nuclear and pure neutron 
matter; the difference between them (dot-dashed lines) 
is the asymmetry energy p7p . By way of comparison, we 
found that the curve for Eq. (|36|) deviates by less than 
2% from that of Eq. ((37|) . which is numerically not signif- 
icant. In all three models, values of rig"*, E^^^, and a^^^ 
are in perfect agreement with the empirical ones [Eqs. ([3]), 
(III), and (O], whereby the standard deviation errors are 
less than 10~^%. All these deviations arise solely from 
the numerical solution of equation ([50]) . In RMF2, the 
saturation point {n^*,E-^^) is reproduced through G^ 
and Gi^ only. 

For isospin symmetric nuclear matter at the saturation, 
bound states of nucleons appear with delicate cancella- 
tion between short-range repulsion and medium-range 
attraction. The system is in its ground state, whereby 
all energy levels are occupied separately by protons and 
neutrons up to the Fermi energy of about 779.3 MeV. 
Outside the saturation region, the system becomes un- 
bound because the binding energy per nucleon becomes 
large (between 300 and 500 MeV) comparing with the 
nucleon mass. 



> 



rmj 



RMF3 (without Gt) 
r; = 0., 0.1, 0.2,...., 1 




r? = 0. (ran — ftp) 



ffh I I I I I I I I n+i 



0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 
Nuclear Density ran [fm^^] 

Figure 4: (Color online) The evolution of the saturation point 
with proton-neutron fraction 77 in RMF3 model. The 7; = 
{ri = 1) solid curve corresponds to isospin symmetric nu- 
clear matter (neutron matter). The rj = 0.1, ...,0.9 dashed 
curves correspond to asymmetric nuclear matter which in- 
creases with rj from bottom to top. We found that nuclear 
matter ceases to saturate at 77 slightly larger than 0.8. 



The evolution of the binding energy per nucleon of 
asymmetric nuclear matter with proton-neutron asym- 
metry r] is illustrated in Fig. [4] only for the nonlinear 
RMF3 model (for reasons that become clear in Sec. lIIip . 
The two solid lines correspond to isospin symmetric and 
pure neutron matter at 77 = and 1, respectively. The 
nine dashed lines correspond to asymmetric nuclear mat- 
ter at 77 = 0.1, ...,0.9. We found that nuclear matter 
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ceases to saturate at rj slightly larger than 0.8, as shown 
in Table ID By way of comparison, this value is slightly 
larger than that obtained by Piekarewicz and Centelles 
using the so-called NL3 and FSUGold model calcula- 
tions |18| . But the conclusion is the same, namely, the 
larger the proton-neutron asymmetry, the smaller the 
saturation density and hence the lower the binding en- 
ergy per nucleon, which means that pure neutron matter 
is an unbound system at all densities. 



etc.), which can be done by fitting to the empirical data. 
In our case, the analytical expression (|40|l includes (a) 
the compression modulus of isospin symmetric nuclear 
matter, which characterizes the small density fluctua- 
tions around the saturation point {n^*,Ef^^), and (b) 
the compression modulus of asymmetric matter, which 
characterizes the slope and curvature of the asymmetry 
energy around the saturation density at given value of rj. 





Results for RMF3 (i.e., without G4) 


P.-N. Asym. 


Saturation Point 


Compress. Mod. 


V 


nr [fm--^] 


£""* [MeV] 


IC [MeV] 


0.0 {Un = Up) 


0.153 


-16.3000 


234.000 (= /Ci,o) 


0.1 


0.151 


-15.9767 


229.343 


0.2 


0.146 


-15.0260 


220.058 


0.3 


0.138 


-13.5048 


205.440 


0.4 


0.127 


-11.5035 


184.452 


0.5 


0.113 


-9.14113 


155.991 


0.6 


0.098 


-6.55846 


128.252 


0.7 


0.081 


-3.91386 


95.8663 


0.8 


0.062 


-1.38341 


59.8410 


0.9 








(?) 


1.0 (np = 0) 








(?) 



Table I: The evolution of the saturation density, the bind- 
ing energy per nucleon and the compression modulus of nu- 
clear matter with proton-neutron asymmetry in the nonlinear 
{a — oj) model, without the quartic self-coupling constant G4. 



In the case of a very asymmetric matter, 0.8 < t] ^ 1, 
the saturation point vanishes systematically, so that we 
cannot obtain a precise value for K, within the {a — uj) 
model (see Table [I]). Nevertheless, we are "curious" to 
know the response of Eq. ^^ for some "hypothetical 
equilibrium" values of tib at r; = 1. This situation is 
quite similar to that of neutron stars in terms of the 
baryon density within an order of magnitude or less. The 
difference, however, lies in the fact that neutron stars are 
bound by gravity, not by the isospin symmetric nuclear 
forces. So, if we consider a superdense neutron matter 
at an "equilibrium density" between 0.8 and 1.0 fm~^, 
expression ([40|) returns very high values of /C = JCpnm as 
presented in Table [III 



We close this section by discussing the behavior 
of the pressure and the compression modulus. For 
isospin symmetric nuclear matter, we obtain from 
Eq. go]) IC^so ^ 556.693 MeV in RMF2, which is much 
larger than the empirical one ([7]). This result is well 
known from the standard Walecka model [ij]. In both 
RMF3 and RMF4 we obtain "exactly what we want," 
ICiso — 234.0 MeV owing (a) to the nonlinear extension 
known originally from Boguta and Bodmer [l7[, (b) to 
the analytic expression for K. given in Eq. (|40|) . and (c) 
to the parametrization of the coupling constants, which 
would be presented in Sec. IIIII 



ns [fm-3] 


m^ [MeV] 


ICpn^ [MeV] 


0.8 


444.724 


8012.41 


0.9 


420.755 


9184.54 


1.0 


399.871 


10342.1 



Table II: The compression modulus of pure neutron matter 
for RMF3 model for some "hypothetical" values of the equi- 
librium density. 



The compression modulus of asymmetric nuclear mat- 
ter decreases with increasing proton-neutron asymmetry 
T] beween 0.0 and 0.8, as shown in the last column of Ta- 
ble H] for RMF3 model. Certain authors like Piekarewicz 
and Centelles [I3 obtained roughly the same results by 
expanding lC{r]) in a power series in 77^, 

IC{r]) = IC^so + ?7^A^sym + 0{lf), 

which is fully analogous to the one for the binding energy 
per nucleon discussed earlier. Other authors like Blaizot 
et al. [8| write /C(r/) empirically as a sum of different types 
of compression modulus (volume, asymmetry. Coulomb, 



Finally, Fig.[S]shows the pressure for isospin symmetric 
nuclear (solid lines) and pure neutron matter (dashed 
lines). Its behavior in RMF3 as well as in RMF4 is almost 
the same owing to the value of the coupling constants 
used here. 
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Nuclear Density riB [fin ^] 

Figure 5; (Color online) The pressure as a function of the 
nuclear density in the RMF approximation to the nonlinear 
(a - oj) models (RMF3 and RMF4) and the linear model 
(RMF2) at zero temperature. The solid lines correspond to 
isospin symmetric nuclear matter and dashed lines to pure 
neutron matter. 



J 



At low density, the nuclear matter is in pure diluted 
phase, similar to a nonrelativistic ideal Fermi gas. In 
this case, expression (f?T|) can be approximated by the 
following fifth-order polynomial of kp 



V 



RMF4 



(UB « n|f*) 



Z k p 



(42) 



At higher density above the saturation, the pressure in- 
creases very strongly and the nuclear medium becomes 
in dense phase (liquid) . In this case, expression ([51]) con- 
verges to 



Pressure/^/' 



iR.MF4 



[uB > n|f*) 




^ . (43) 



Pressure/^ /^/^ 



PrcssurcXj 



correction 1st. order /^ 



negligible 



Here the leading term (which is the first) comes from the 
strong repulsive forces caused by the vector a; meson ex- 
change. Afterward, small negative terms follow owing to 
attractive forces caused by the scalar a meson exchange 
(there is always exchange even at higher density) , and the 
last three positive terms are their first relativistic correc- 
tions, which slightly increase the strength of the scalar 
interactions. Therefore, the second and third terms give 
the difference between the linear and the two nonlinear 
(ct — cj) models. 



III. ANALYTIC PARAMETRIZATION OF THE 
(a - a;) MODEL 



First, let us collect the five empirical quantities for 
isospin symmetric nuclear matter discussed in the Intro- 



duction 

the saturation density: n|°* = 0.153 fm"^, (44a) 

the binding energy per nucleon: E^^* = —16.3 MeV, 

(44b) 

the asymmetry energy coefficient: afy^ = 32.5 MeV, 

(44c) 

the nucleon effective mass: 0.7 < m^/m-^ < 0.8, 

(44d) 

the compression modulus: 230.0 < KLiso S 250.0 MeV. 

(44e) 

In this section, we answer the following questions 

(i) Which values of the coupling constants of the 
(a — cj) model do we have to input in the EoS in or- 
der to reproduce the five empirical quantities (|44p ? 



(ii) Is it possible to reproduce them without a four-body 
self- interaction? 

The input parameters for us here will be the five empirical 
quantities (|44|) . while the output parameters will be 

• G„, G^, Gp, G3, and Gi in RMF4, 
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• Ga;Gu;,Gp, and G3 in RMF3, and 

• Ga, G^, and Gp in the RMF2 model. 

Let us rewrite all EoS for cold isospin symmetric nuclear 
matter; after evaluating the integrals in Eqs. ([50)1 . ([5T|) . 



and (|33p . we have the following. 



J 



The nucleon effective mass 
G. 



m-^ = Ton 
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where ICiso is deduced from (j40p for ?] = 0. Further we have 



"B = -^^Fi 



E- = 



B,. 



A 



cRMFi 



Ton, ctq = Ton - Ton 
I 



and Ep = Hq = y k 



"F 



In the next three sections, all quantities ne, Eiso, kp, 
m^, Ep, and asym are evaluated at the saturation (ne = 

nff*, £;,,„ = £;-*, etc.). 



A. Parametrization of the RMF4 Model 

The isospin coupling constant can be immediately ob- 
tained from Eq. (|48|) as 



Gp - ^^ ^asym g^. 



(50) 



By adding the EoS (|46cp to Eq. (|47cp and requiring that 
V = MeV/fm'^, we obtain the repulsive coupling con- 
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stant 



tx^ — 



Eiso + rrif^ — Ep 



riB 



(51) 



Now, to determine G^, G3, and G^, wc need three equa- 
tions, for example, Eqs. (gS]), dH]), and ((i7c)l [or pSc)) ]. 
From Eqs. (|45|) and (|47cp . wc obtain 



ce3 ce4 
G3 G4 



ce5. 



(52) 



Further, from Eqs. (|451) . and p9|) wc obtain 
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Finally, from Eqs. ([5^ and ([55)1 we deduce G^, G3 and 
G4 in terms of the empirical quantities tib , Eiso i "t-n i ^^'^ 
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* \i 



kpim^ - Ton) 



^ {kl + 3m;f) ( 6E*pkl + tt^ubIC^so - 6k%{E,so + ton) ) - 3 (^fc^^^) 



o XT'* 2 *2 
SA^TOnTIn 



7r2(3fc^ - E*pIC,so) E*pk%{E*p - E,, 

Ub 



Ton 



In 



E*p + kF 



+ 



(54) 



G3 = <^7r2£:|,(TON-TOj:j)3 



n\3kj. - £:>/C,,o) E*pk%{E*p - ^,,0 - ton) 



15Ep kp{mj<i — m^)ml 



TT^ {3k% - EpICiso) ~ Epkp{E*p - Eiso - ^n) 



ne 



3/2 /^(ton — TOn)^ 



3 {nkpEpf + ^> (fc|, + 3toJ;,2) f^k%{E,so + TON - £^1^) - ^r^JC^soj 



2E*p 



6 

nB 



—E*pk],{E*p - £;,,„ - ton) - TT^ (3fc| - £:|,/C,.o) 



QTT^UBiEiso + to.n) + -B|, 2k% - 3k pm*^ - Qn^UB 



3m^E*pm:^ {3mT^ - m^) 



—Efkl - —E*pk%{E,so + ton) - n^ (3kl - E*plC,s 
riB riB 



In 



E*+k, 



(55) 



14 



< 8E'pkp{mj<s — m^)m 



n^ (3k% - E*pJCiso) EpkpiEp - Eiso - m^) 



UB. 



-Epkp(E*p - Eiso - ?«n) - S-n-'kp + n-EpJCiso 



2kp{mj^ — 771^) 
E 



3 {TTkpEpy + E*p (k'p + 3m;f) —ep{E,so -E*p + mN) - rr'lC^so 



^2 (3fc^ - E*pJC.,so) E*pkl{E*p - E,so - mN) 



67r 7iB(£'iso + "^n) + Ep 2k p ~ 3kpm'^ - 67r tib 



Epm^ ( Smjvj — 4ton"In + m* 



—E*p^kl - —E*pkl{E,so + mN) - 3^2fc2^ + 7r^E*pIC^so 
riB nB 



In 



S*^ + fe, 



.(56) 



As we see, they are expressed as rational functions, 
{...}/{...}. The relation between G^ and the nonlinear 
coupling constants is given by 



G. = 



c^o 



al/Gs - a^o/Gi 



(57) 



The results are shown in Table IIIII for different values 
of ICiso between 230 and 250 MeV and m'^/m-^ between 
0.7 and 0.8. In each line, the set of coupling constants 
yields the saturation density (|44a|) . the binding energy 
per nucleon at saturation (j44bp and the asymmetry en- 
ergy cocfhcicnt (|44cp with standard deviation errors less 
than 10~^% arising from the numerical solution of equa- 
tion Eq. ([45]) only. Thus, G3 and G4 together provide 
high accuracy in determining of the five empirical quan- 
tities (144]). This accuracy was also confirmed by other 
authors [19l |20[ by investigating the ground state prop- 
erties of finite nuclei through accurately calibrating pro- 
cedures, such as the NL3 and FSUGold models. 
However, for 230.0 < Jd^o ^ 232.7 MeV, we found 
that the quartic self-coupling constant G4 is negative at 
all values of rri^/m^. Further, at Kiso = 234.0 MeV 
(the value that we have supported in this work [^) 
and for m^/mN ^ 0.79, G4 is always negative, while 
for mJj/mN > 0.79 it takes very large positive values 
(above 1987.3). Notice that G4 should be less than 
1000 to have an effect of roughly 22.1% on the bind- 
ing energy per nucleon at the saturation. Here we found 
that it is only possible to get G4 between and 1000 if 
0.731 < m;j/77iN < 0.8 and 235.3 < JCiso < 250.0 MeV. 
Of course, mathematically it does not matter if G4 is neg- 
ative or positive, because all sets of coupling constants 
in Table Uni which fulfill Eqs. dSO]), ([531) and dSZ]) at sat- 
uration, yield the same behavior of the EoS at all nu- 
clear densities (sec Fig.s [21 [21 and [5] for RMF4, where 
G4 = -254.991). 





Nr. 


rriN 


[McV] 


(GcV-2j 


G^ 

[GcV-2] 


[GcV-2] 


G3 

[GcV-l] 


G4 


1 


0.71 


230.0 


311.430 


176.413 


114.672 


204.403 


-182.178 


2 


0.71 


232.0 


310.946 


176.413 


114.672 


206.731 


-186.768 


3 


0.72 


230.0 


303.947 


168.962 


115.960 


188.641 


-174.967 


4 


0.72 


232.0 


303.442 


168.962 


115.960 


190.892 


-179.960 


5 


0.72 


234.0 


302.937 


168.962 


115.960 


193.204 


-185.258 


6 


0.78 


234.0 


257.169 


123.993 


123.119 


116.736 


-254.991 


7 


0.73 


230.0 


296.465 


161.498 


117.220 


173.707 


-169.047 


8 


0.73 


232.0 


295.936 


161.498 


117.220 


175.887 


-174.583 


9 


0.73 


234.0 


295.408 


161.498 


117.220 


178.126 


-180.508 


10 


0.73 


236.0 


294.880 


161.498 


117.220 


180.429 


-186.863 


11 


0.74 


230.0 


288.972 


154.021 


118.452 


159.642 


-164.922 


12 


0.74 


232.0 


288.417 


154.021 


118.452 


161.755 


-171.233 


13 


0.74 


234.0 


287.863 


154.021 


118.452 


163.930 


-178.063 


14 


0.74 


236.0 


287.310 


154.021 


118.452 


166.170 


-185.479 


15 


0.75 


230.0 


281.456 


146.532 


119.657 


146.478 


-163.434 


16 


0.75 


232.0 


280.872 


146.532 


119.657 


148.533 


-170.930 


17 


0.75 


234.0 


280.289 


146.532 


119.657 


150.653 


-179.168 


18 


0.75 


236.0 


279.706 


146.532 


119.657 


152.839 


-188.264 


19 


0.76 


230.0 


273.900 


139.030 


120.836 


134.247 


-166.160 


20 


0.76 


232.0 


273.283 


139.030 


120.836 


136.256 


-175.643 


21 


0.76 


234.0 


272.666 


139.030 


120.836 


138.332 


-186.304 


22 


0.77 


230.0 


266.280 


131.517 


121.990 


122.983 


-176.534 


23 


0.77 


232.0 


265.625 


131.517 


121.990 


124.961 


-189.847 


24 


0.79 


246.0 


244.787 


116.459 


124.225 


122.098 


682.043 


25 


0.79 


248.0 


244.056 


116.459 


124.225 


124.937 


477.204 


26 


0.79 


250.0 


243.327 


116.459 


124.225 


127.922 


366.713 


27 


0.80 


236.0 


240.243 


108.914 


125.307 


101.935 


777.923 


28 


0.80 


238.0 


239.447 


108.914 


125.307 


104.286 


482.989 


29 


0.80 


240.0 


238.654 


108.914 


125.307 


106.756 


349.882 


30 


0.80 


242.0 


237.863 


108.914 


125.307 


109.355 


274.087 


31 


0.80 


244.0 


237.074 


108.914 


125.307 


112.093 


225.146 


32 


0.80 


246.0 


236.288 


108.914 


125.307 


114.982 


190.935 


33 


0.80 


248.0 


235.505 


108.914 


125.307 


118.035 


165.674 


34 


0.80 


250.0 


234.724 


108.914 


125.307 


121.266 


146.257 



Table TIL Coupling constants of the nonlinear (ct — ui) model 
for nuclear matter, RMF4, for different values of the nucleon 
effective mass and the compression modulus, that reproduce 
the saturation density, the binding energy per nucleon at sat- 
uration and the asymmetry energy coefficient. 
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Although, negative values of G4 do not constitute good 
news for the nonlinear (a — w) model, we should not dra- 
matize this for the reason that the scalar a field is treated 
as a classical mean field, which is not realistic. Further, 
the (tr — uj) model treats baryons and mesons as elemen- 
tary degrees of freedom; this can be one of the reasons 
why various values of G4 are negative at given values of 
K-iso and m^. To be closer to realistic descriptions, the 
coupling constants should, for example, depend on the 
momentum transfer between nucleons. In other words, 
negative coupling constants will be a serious problem 
when (a) we examine the {a — tu) model with momentum- 
dependent coupling constants of the meson fields, that 
is, by incorporating form factors into the nucleon-meson 
couplings, or (b) we derive the EoS for nuclear matter 
directly from QCD, which is, however, not achieved so 
far. 



B. Parametrization of the RMF3 Model 



In RMF3 model we take G4 — > ±cx3. To give a pos- 
itive answer to the second question of this section, the 
coupling constant G3 must 



Nr. 


JTlN 


[MeV] 


G. 

[GeV-2] 


[GeV-2] 


G, 

[GeV-2] 


G3 

[GeV-^l 


1 


0.775 


253.862 


254.463 


127.757 


122.558 


147.074 


2 


0.780 


249.712 


251.748 


123.993 


123.119 


136.037 


3 


0.781 


248.878 


251.212 


123.240 


123.231 


133.913 


4 


0.782 


248.043 


250.679 


122.487 


123.342 


131.815 


5 


0.784 


246.368 


249.619 


120.981 


123.564 


127.699 


6 


0.786 


244.687 


248.569 


119.474 


123.785 


123.687 


7 


0.787 


243.844 


248.048 


118.720 


123.896 


121.718 


8 


0.788 


242.999 


247.530 


117.967 


124.005 


119.775 


9 


0.790 


241.304 


246.502 


116.459 


124.225 


115.963 


10 


0.793 


238.746 


244.983 


114.196 


124.552 


110.426 


11 


0.794 


237.889 


244.483 


113.442 


124.660 


108.627 


12 


0.795 


237.030 


243.986 


112.687 


124.769 


106.851 


13 


0.797 


235.305 


243.003 


111.178 


124.985 


103.368 


14 


0.798 


234.450 


242.521 


110.424 


125.092 


101.656 


15 


0.798 


234.001 


242.276 


110.046 


125.146 


100.815 


16 


0.799 


233.570 


242.033 


109.669 


125.200 


99.975 


17 


0.800 


232.710 


241.559 


108.914 


125.307 


98.307 



Table IV: Coupling constants of the nonlinear (a — lo) model, 
RMF3, for different values of the nucleon effective mass and 
the compression modulus, that reproduce the saturation den- 
sity, the binding energy per nucleon at saturation and the 
asymmetry energy coefficient. 



52]) and (1531 ■ which is 



• verify simultaneously Eqs. 
the best condition, or 



• verify Eq. (|52p and "approximately" Eq. ([53]), or 



verify Eq. (|53p and "approximately" Eq. ([5^ . 



If 6*3 =: ce3/ce5, then mjj and Eiso will be reproduced, 

but there is no guarantee for K-iso- If G3 = cq/cIG, then 
m^ and ICiso will be reproduced, but there is no guaran- 
tee for Eiso- Thus, let us define the numerical parameter 



diffe 



def 



cl6 



ce3 

ce5 



(58) 



So, if diffe — > and G3 verify Eq. (|52]) or (|53l), then 
the five empirical quantities (|44l) can well be reproduced 
in RMF3 model. The smaller the parameter diffe is, 
the more precisely we obtain our desired results. In Ta- 
ble |lVl we present some results for diffe ^ 10^^ and by 

requiring that G3 = ce3/ce5. 



Certainly we can obtain roughly the same sets of cou- 
pling constants by using Eqs. ([45]) . (|49|) . and (|46cp [in- 
stead of Eq. (|47c|) ]. and then we can parametrize the 
RMF3 model by defining the following numerical param- 
eter: 



dif f p "^^ 


(To CpZ 

cl6 cp5 


= —cp2 , 


ctq 
3 ' 


= ^Gunl+ 


'2 
-E*pk%, - E*pkFm*^ - m*^ In 





(59) 



El 



Pf 



where 
cp5 : 
cp3 : 

cp2: 

1 

4^ 



Finally, if we compare the results of Table Hill with those 
of Table IIVI we see that the accuracy in terms of the 
nucleon effective mass and/or the compression modulus 
obtained earlier in RMF4 becomes more difficult in the 
RMF3 model. This is the price we have to pay in RMF3, 
which is of minor importance as long as the empirical 
values of m^ and IC cannot be estimated (at present) to 
an accuracy of better than ^ 0.02ton and '^ 0.2 MeV, 
respectively. 
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C. Parametrization of the RMF2 Model 

To be in the linear (a — uj) model, we can take G3 — > 
+CXD and G4 — > ±cx). The scalar coupling constant be- 
comes Ga = (To/'^cr- This, together with G^j and Gp, 
yields just riB, Eiso and agym at saturation. Some values 
are listed in Table IVl 



Nr. 


?7i]\j/mN 


;C,.o[MeV] 


Ga[GeV-^] 


G^[GeV-'] 


Gp[GeV-'] 


1 


0.536 


576.368 


398.114 


303.302 


86.608 


2 


0.537 


570.700 


397.293 


302.592 


86.806 


3 


0.538 


565.052 


396.472 


301.881 


87.004 


4 


0.539 


559.441 


395.651 


301.171 


87.201 


5 


0.540 


556.693 


394.831 


300.460 


87.398 


6 


0.541 


548.265 


394.011 


299.748 


87.594 



Table V: Coupling constants of the linear {a — uj) model, 
RMF2, that reproduce the saturation density, the binding 
energy per nucleon at saturation and the asymmetry energy 
coefBcient, but not the compression modulus and the nucleon 
effective mass. 

Thus, we clearly confirm here the results obtained by 
the standard Walecka model M. 



IV. CONCLUSIONS 

From all these results, we can conclude that it is 
helpful [i^ [20], but not necessary, to introduce the 
four-body self-interaction into the linear {a — cu) model, 
as demonstrated in the present work. As long as the 
empirical values of the nucleon effective mass and the 
compression modulus cannot be estimated with high ac- 
curacy, a three-body self-interaction would be sufficient 
(and of course necessary) , so that we can prevent getting 
a negative coupling constant. Briefly, G4 almost always 
plays the role of a negative fine-tuning parameter 
of the nonlinear (c — uj) model for bulk nuclear matter. 

Furthermore, we expect that the same conclusions 
can be reached for finite nuclei, because the Lagrangian 
densities are not too different (see, for example, Ref |19|'): 
that is, the only difference is that we have to include 
the coupling of protons to the electromagnetic field 
and a derivative surface term. We point out that the 
parametrization of the (a — u) model was also made by 
Glendenning [16| for isospin symmetric nuclear matter 



(t] ~ 0) and also found that the quartic self-coupling 
constant is negative. It is now, in principle, possible to 
determine the coupling constants of the {a — uj) model 
analytically for asymmetric nuclear matter (0 < 77 < 1), 
because we have now an analytic expression for the 
compression modulus given in Eq. (|^0)) . This calculation 
makes the effective model especially well defined in 
terms of the empirical properties and the effect of the 
proton-neutron asymmetry on the values of the coupling 
constants. This study would be most welcome. 

In summary, we have discussed the empirical proper- 
ties of bulk nuclear matter obtained from modern scatter- 
ing experiments. We have calculated the EoS for infinite 
nuclear matter within the linear and nonlinear (a — uj) 
model at finite temperature. Afterward, our attention is 
focused on the zero-temperature EoS. Within the non- 
linear model, we have presented an analytical expression 
for the compression modulus and investigated the EoS at 
all nuclear densities and different proton-neutron asym- 
metry 77. We have found that nuclear matter ceases to 
saturate at 77 slightly larger than 0.8. 
In Sec. mil we have developed an analytical method to 
determine the strong coupling constants from the EoS 
for isospin symmetric nuclear matter, which allow us to 
reproduce all the saturation properties. This method has 
allowed us to interpret the quartic self-coupling constant 
as a negative, or positive and very large fine-tuning pa- 
rameter of the nonlinear {a — to) model. We have clearly 
confirmed |14{ that the linear model cannot reproduce 
the nucleon effective mass and the compression modulus, 
so a nonlinear extension introducing self-coupling terms 
of the scalar meson field seems to be necessary. However, 
owing to this analytical method, we have shown that it 
is possible (a) to investigate the EoS in terms of ne and 
1] and (b) to reproduce all the five saturation properties 
without the quartic self-coupling constant. We have thus 
concluded that the latter is not necessary in the {a — uj) 
model. 
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